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' Abstract 

The purpose of this paper is to investigate the interplay arising between max algebra, convexity and 
^\ , scaling problems. The latter, which have been studied in nonnegative matrix theory, are strongly 

04 ' related to max algebra. One problem is that of strict visualization scaling, defined as, for a given 

nonnegative matrix A, a diagonal matrix X such that all elements of X~^AX are less than or equal 
to the maximum cycle geometric mean of A, with strict inequality for the entries which do not lie on 
critical cycles. In this paper such scalings are described by means of the max algebraic subeigenvectors 
and Kleene stars of nonnegative matrices as well as by some concepts of convex geometry. 
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^ ■ 1 Introduction 

(N 

0\ ■ 

0\ ■ The purpose of this paper is to investigate the interplay arising between 

^ max algebra, convexity and matrix scaling. A nonnegative matrix A is 

O ' called visualized if all its elements are less than or equal to the maximum 

00 ' 

O ■ cycle geometric mean X{A) of A, and it is called strictly visualized if, fur- 

ther, there is strict inequality for the entries which do not lie on critical 
^ ■ cycles. Given a nonnegative matrix A, the chief aim of this paper is to 

5^ ; identify and characterize in several ways diagonal matrices X with a pos- 

itive diagonal for which X^^AX is strictly visualized, see Theorems 3.3, 
3.7, 4.2 and 4.4. 

In Section 2, we revisit and appropriately summarize the theory of max 
algebraic eigenvectors and subeigenvectors, and some properties of Kleene 
stars. 
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Sections 3 and 4 contain our principal results. In Section 4 our chief 
tool is the Kleene star ^4* of A (defined for a definite matrix), and the max 
algebraic cone V*{A). The latter consists of the subeigenvectors of A for 
the eigenvalue X{A) or, equivalently, of the eigenvectors of A*. We call 
V*{A) the subeigencone of A. It is also a convex cone. Diagonal matrices 
X corresponding to vectors x in its relative interior of the subeigencone 
are precisely the matrices X that strictly visualize A, see Theorem 3.7. 
Among those vectors x are all linear combinations of the columns of A* 
with positive coefficients, see Theorem 3.3. 

While in Section 3 our approach is convex geometric, the main idea of 
Section 4 is to start with a strictly visualized matrix and to describe all 
strict visualizers in matrix theoretic terms, see Theorem 4.2. We also show 
that the dimension of the linear hull of the subeigencone V*{A) equals 
the number of components of the critical graph of the Kleene star ^4*, see 
Theorem 4.4. At the end of the section we show by example that the max 
algebraic dimension of V*{A) may exceed its linear algebraic dimension. 

The interplay between max algebra (essentially equivalent to tropical 
algebra) and convexity, here explored via visualization, is also important 
for tropical convexity, see the papers [15], EH, [32], among many others. We 
also note that visualization scalings can be important for max algebra, due 
to the connections with the theory of — 1 matrices that they provide. See 
[m, dZI, [39] for recent developments and apphcations of this idea. 

2 Eigenvectors and subeigenvectors 

By max algebra we understand the analogue of linear algebra developed 
over the max-times semiring Mmax,x which is the set of nonnegative numbers 
M+ equipped with the operations of "addition" a © 6 := max(a, b) and the 
ordinary multiplication a b := a x b. The operations of the semiring are 
extended to the nonnegative matrices and vectors in the same way as in 
conventional linear algebra. That is if ^4 = (a^j), B = [bij) and C = (c^) 
are matrices of compatible sizes with entries from R+, we write C = A(B B 
if Cij = ttij^bij for aU ij and C = A0B if Cij = J2k O'lkbkj = uidiXkiaikbkj) 
for alH, j. If Of G M+ then a A = {aaij). We assume everywhere in this 
paper that n > 1 is an integer. P„ will stand for the set of permutations of 
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the set {1, n}, and the sets hke {1, . . . , m} or {1, ... , n} will be denoted 
by [m] or [n] , respectively. If ^4 is an n x n matrix then the iterated product 
A<S) A<S) ... 'S) A in which the symbol A appears k times will be denoted by 
AK 

Max algebra is often presented in settings which seem to be different from 
IRmax,x, namely, over the max-plus semiring Mmax,+ = (M U {— (X)},© = 
max, = +) and the min-plus (or tropical) semiring Mmin,+ = (M U 
{+oo},© = min, = +). The semirings are isomorphic to each other 
and to Mmax,x- In particular, x i-^ exp(a:) yields an isomorphism between 

IK.max,+ Q-Ild M-niax, x • 

Let A = (aij) G M"^". The max algebraic eigenproblem consists in 
finding A G M+ and x G M'|\{0} such that A0 x = Xx. If this equation is 
satisfied, then A is called a max algebraic eigenvalue of A and x is called a 
max algebraic eigenvector of A associated with the eigenvalue A. 

We will also be interested in the max algebraic subeigenvectors associated 
with A, that is, x G such that A^ x < Xx. Their first appearance in 
max algebra seems to be [22] Ch. IV and [23]. For a more recent reference, 
see generalization of the max-plus spectral theory fl\ , where they are called 
super-eigenvectors. 

Next we explain two notions important for both the eigenproblem and 
the subeigenproblem: that of the maximum cycle mean and that of the 
Kleene star. 

Let A = {a^j) G M^^". The weighted digraph Da = {N{A), E{A)), with 
the set of nodes N{A) = [n] and the set of edges E{A) = N{A) x N{A) 
with weights w{i^j) = a^j, is called the digraph associated with A. Suppose 
that TT = (^1, Zp) is a path in Da, then the weight of tt is defined to be 
w{7r, A) = tti^i^tti^i^ . . . ai^_^i^ if p > 1, and 1 if j9 = 1. If = ip then tt is 
called a cycle. A path tt is called positive if w{7t,A) > 0. A path which 
begins at i and ends at j will be called an i ^ j path. The maximum cycle 
geometric mean of A, further denoted by X{A), is defined by the formula 

X{A) = max jj,{a^ A), 

a 

where the maximization is taken over all cycles in the digraph and 

li{a,A) = w{a,Ay/'' 
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denotes the geometric mean of the cycle a = (ii, 

If the series / ^4 © ^4^ © . . . converges to a finite matrix, then this 
matrix is called the Kleene star oi A and denoted by A* = {alj)- The next 
proposition gives a necessary and sufficient condition for a matrix to be a 
Kleene star. 

Proposition 2.1 ]J] Let A = (a^j) G W;^^ . The following are equivalent: 

1. A is a Kleene star; 

2. A' = A; 

3. A'^ = A and an = 1 for all i = 1, . . . ,n. 

The next theorem explains some of the interplay between the maximum 
cycle geometric mean X{A), the Kleene star A*, and the max algebraic 
eigenproblem. 

Theorem 2.2 fl, d 23, gg/ i^ei ^ G M^^". Then 

1. the series / © ^4 © ^4^ © . . . converges to a finite matrix A* if and only 
if A(^) < 1, and then A* = I ® A ® A^ ® . . . ® A^-^ and A(^*) = 1; 

2. is the greatest max algebraic eigenvalue of A. 

This theorem shows great similarity between max algebra and nonnega- 
tive linear algebra. However, it also reveals a crucial difference: the series 

. converges also if A(74) = 1. 

A G R"^*^ is called irreducible if for any nodes i and j in Da a positive 
i ^ 3 path exists. 

Proposition 2.3 ^ [13/ If A is irreducible and < 1, then A* has all 
entries positive. 

More generally, it is important that Kleene stars accumulate the paths 
with greatest weights. Namely, if i j then a*j = maxi(;(7r, ^4) where tt 
ranges over paths from i to j. 

Matrices with X{A) = 1 are called definite. 

Results involving a Kleene star A* will be stated for definite matrices. 
There is no real loss of generality here in the case of matrices A with 
X{A) > 0. Indeed, for any such A we have that X{aA) = aX{A), and if 
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a > 0, then any eigenvector of A associated with X{A) is also an eigenvector 
of aA associated with X{aA) and conversely. Hence if X{A) > 0, then the 
eigenproblems for A and A/X{A), which is definite, are equivalent. 

Note that X{A) = implies that A contains a zero column, and then 
eigenvectors and subeigenvectors are just vectors x satisfying Xi = when- 
ever the corresponding column A.i ^ 0. In what follows, we will not treat 
this trivial case and we will always assume that X{^A) > 0. 

The spaces that we consider in max algebra are subsets of closed 
under componentwise maximization ©, and scalar multiplication. They are 
called max cones, due to the apparent analogy and important connections 
with conventionally convex cones in MJ^. 

The set of subeigenvectors of A associated with X{A) will be denoted 
by V*{A). The set of eigenvectors associated with X{A) will be denoted 
by ^(^4). Both sets are max cones, and hence ^(^4) will be called the 
eigencone of A, and V'^{A) will be called the subeigencone of A. Next we 
study some simple relations between ^(^4) and V*{A). The first one is 
immediate. 

Proposition 2.4 V{A) C V*{A). 

Further we denote by span^(74) the max algebraic column span of A, 
which is the set of max combinations {Y^f XiA.i, Xi G M+} of the columns 
of A. Note that V{A) C sp8in^{A) for any matrix A. 

Proposition 2.5 // A is definite, then V'{A) = V{A') = V*{A*) = 
span^{A*). 

Proof. First note that by Theorem 12.21 if X{A) = 1 then A* exists and 
A(^*) = 1. Now we show that V*{A) = V{A*). Suppose that A* ®x = x, 
then A <S) X < X, because A < ^4*. If A<S) x < x, then {I ® A) <S) x = x and 
also A* <S) X = X, since A"^ x < x for any m (due to the monotonicity 
of matrix multiplication). As (^4*)* = A* by Prop. 12.11 we also have that 
= V{A'). 

We show that V*{A) = spangg(^*). As A^A* < A*, each column of A* is 
a subeigenvector of A, hence span^(74*) C V*{A). The converse inclusion 
follows from V*{A) = V{A*) and the inclusion V{A*) C span^(74*). ■ 
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A matrix A will be called strongly definite, if it is definite and if all its 
diagonal entries equal 1. Note that any Kleene star is strongly definite by 
Prop. O 

Proposition 2.6 For A a strongly definite matrix, V{A) = V*{A). 

Proof. To establish V{A) = V*{A), it is enough to show ^{A) C V{A), 
as the converse inclusion is trivially true. Take y G V*{A). We have that 
aijUj^Ut < Hi which is equivalent to aijUj^Ut = Hi, so y e V{A). 

■ 

By the above propositions, the subeigenvectors of A, and in the strongly 
definite case also the eigenvectors of A, are described as the vectors from 
the max algebraic column span of A*, which we call Kleene cone. 

More generally, a set S is called a generating set for a max cone 
written K = span^(5'), if every vector y E K can be expressed as a max 
combination y = A^x* of some elements G 5, with Xi > 

for i G [m]. A set S is called a (weak) basis for K if span^(5') = K and 
none of the vectors in S can be expressed as a max combination of the 
other vectors in S. A vector y E K is called a max extremal of if 
y = u (B u,w E K implies that y = u or y = w. The set of max 
extremals u of K scaled with respect to the max norm, which means that 
||m|| = maxjWj = 1, will be denoted by ext©(i^'). We have the following 
general result describing max extremals of closed max cones. 

Theorem 2.7 /771 \^ If K C is a closed max cone, then the set 
extQ){K) is non-empty and it is the unique scaled basis for K. 

li K = span^(74) for some matrix A, then K is closed, so the set 
ext®(span(^(74)) denoted by ext®(74) for brevity, is non-empty and con- 
stitutes the unique scaled basis for span^(74). In this case the vectors of 
QyiiQ){A) are some of the columns of A scaled with respect to the max norm. 

Next we describe the eigencone and the subeigencone of ^4 G R*]:^", and 
the sets of their scaled max extremals, in the case \{A) > 0. For this we 
will need the following notions and notation. The cycles with the cycle ge- 
ometric mean equal to \{A) are called critical, and the nodes and the edges 
of Da that belong to critical cycles are called critical. The set of critical 
nodes is denoted by Nc{A), the set of critical edges is denoted by Ec{A), 
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and the critical digraph of A, further denoted by C{A) = {Nc{A) , Ec{A)) , 
is the digraph which consists of all critical nodes and critical edges of Da- 
All cycles of C (A) are critical |4j . The set of nodes that are not critical is 
denoted by Nc{A). By C*{A) we denote the digraph with the set of nodes 
[n] and the set of edges E*{A) containing all the loops {i, i) for i G [n] and 
such that {i^i) G El{A)^ for i ^ j, if and only if there exists an i ^ j path 
(zi, . . . , -ip) in C{A). The following theorem describes both subeigencone 
and eigencone in the case when A is definite. For two vectors x and we 
write X y \i X = \y ioT \ > 0. 

Theorem 2.8 Let A G W]^^ he a definite matrix, and let M{A) denote a 
fixed set of indices such that for each strongly connected component ofC{A) 
there is a unique index of that component in M{A). Then A* is strongly 
definite, and 

1. the following are equivalent: {i,j) G Ec{A), aija*^ = a*^ for all k G [n], 
alj = al^Qij for all k e [n]. 

2. the following are equivalent: {i,j) G E*{A), A*- ~ A*j, A*, ~ A*,; 

3. any column of A* is a max extremal o/ span^(74*); 
4- V{A) is described by 



and e'x.i,^{V {A)) is the set of scaled columns of A* whose indices belong 
to M{A); 

5. for any y G V*{A) and any {i^j) G Ec{A) we have aijyj = yi; 

6. V*{A) is described by 



and e'xX,Q{V* {A)) = ext©(74*) is the set of scaled columns of A* whose 
indices belong to M{A) U Nc{A). 



V{A) 
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Proof. Statements 1.-4. are well-known P, [13, [H [221 



We show 5.: By Prop. 12.51 any y G V*{A) is a max combination of the 
columns of A*. Let j) G Ec{A)^ then part 1. implies that aijZj = Zi for 
any z = A\^ /c G [n]. As is a max combination of all these, it follows 
that aijUj = y^. 

We show 6.: By Prop. 12.51 we have V*{A) = span^(^*) and any column 
of A* is a max extremal of spdLnQ{A*) by part 3. By 2. we have that 
A*- ~ A*j if and only if G E*{A), hence aU the columns in M{A) are 
independent max extremals and any other columns with indices in Nc{A) 
are proportional to them. Also note that there are no edges G E*{A) 
such that i ^ Nc{A) or j ^ Nc{A) except for the loops, and therefore all 
columns in Nc{A) are also independent max extremals. ■ 

The number of connected components of C{A) will be denoted by n{C{A)). 
For a finitely generated max cone K the cardinality of its unique scaled 
basis will be called the max algebraic dimension of K. Parts 4. and 6. of 
Theorem 12.81 yield the following corollary. 

Proposition 2.9 For any matrix A G M"^'^ with X{A) > we have that 
the max algebraic dimension ofV{A) is equal to n{C{A)), and the max 
algebraic dimension ofV*{A) is equal to n{C{A)) + |A^c(^)|- 

For X G denote by diag(x) the diagonal matrix with entries 6ijXi, 
for i^j G [n], where 6ij is the Kronecker symbol (that is, 6ij = 1 if i = j 
and 6ij = if z 7^ j). Note that the max algebraic multiplication by a 
diagonal matrix is not different from the conventional multiplication, and 
therefore the notation will be omitted in this case. If x is positive, 
then X = diag(a:) is invertible both in max algebra and in the ordinary 
linear algebra, and the inverse has entries 6ijX^^, for i,j G [n]. The 
spectral properties of a matrix A do not change significantly if we apply 
a diagonal similarity scaling A X~^AX, where X = diag(x), with a 
positive X G M" . 

The following proposition follows very easily from results in the diagonal 
scaling literature, see e.g. Remark 2.9 of [IB] 



Proposition 2.10 Let A G Rl'''' and let B = X^AX, where X = 
diag(a:), with positive x G M'l. Then 
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1. w{a,A) = w{a,B) for every cycle a, hence X{A) = X{B) and C{A) = 
C{B); 

2. V{A) = {Xy I y e V{B)} and V*{A) = {Xy \ y e V*{B)} 

3. A is definite if and only if B is definite, and in this case B* = X~^A*X. 

3 Subeigenvectors, visualization and convexity 

We call X G Wl a nonnegative linear combination (resp. a log-convex 
combination) of y^, . . . ,y'^ E , if x = Xiy^ with > (resp. x = 
UZiiy')^' with A, > and J2Zi ^^ = 1. and both power and multiplication 
taken componentwise). The combinations are called positive if Aj > for 
all i. A set ii' C is called a convex cone (resp. a log-convex set), if it is 
stable under linear combinations (resp. under log-convex combinations). 

In max arithmetics, aQ)b < c is equivalent to a < c and b < c. Using this, 
one can write out a system of very special homogeneous linear inequalities 
which define the subeigencone of A, and hence this cone is also a convex 
cone and a log- convex set. 

Proposition 3.1 Let A e R^'''' and X{A) > 0. Then V*{A) is a max 
cone, a convex cone and a log- convex set. 

Proof. We have that 

V*{A) = {y\A0y< XiA)y} = {y \ J]®a,,i/, < XiA)y, Vz} = 

j 

= {y I (^ijyj < X{A)y, Vi, j}. 

Each set {y \ aijyj < X{A)yi} is a max cone, a convex cone and a log-convex 
set, hence the same is true about V*{A), which is the intersection of these 
sets. ■ 

The log-convexity in (]R+\{0})" (i.e. in the max-times setting) corre- 
sponds to the conventional convexity in R"^ (i.e., the max-plus setting or 
the min-plus setting). We also note that {y \ aijyj < X{A)yi} and hence 
V*{A) are closed under some other operations. In particular, V*{A) is 
closed under componentwise p-norms defined by {y®pz)i = {y^ + z^)^/^ 
for p > 0. 
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Prop. 13. II raises a question whether or not there exist max cones contain- 
ing positive vectors, which are finitely generated and convex, other than 
Kleene cones. The results of [32] suggest that the answer is negative. 

Let K he di convex cone, then y E K is called an extremal of K if and 
only if y = Xu + fiv, where u^v G implies y ^ u (and hence also y v). 
The set of scaled extremals of K will be denoted by ext(i^'). 

Proposition 3.2 Let A e M!^^" and \{A) > 0, then ext©(y*(^)) C 
ext(y*(^)). 

Proof. Without loss of generality we assume that A is definite. By The- 
orem 12.81 part 6., extQ){V* (A)) is the set of scaled columns of A*, after 
eliminating the repetitions. As a*f^ali^ = a*^, for all i^k G [n], we have 
that the x := A\ satisfies a*^Xk = Xi for all i e [n]. As V{A) = V*{A*) 
by Proposition [231 we have that a*i^Zk < Zi for any z G ^{A) and all 
i G [n], implying that if x = Xz^ + fiz'^ with z\ G V*{A), then a'^f^z^ = zf 
for all i G [n] and s = 1,2. Hence 2^ ~ a: and ~ x meaning that 
X G ext(y*(^)). ■ 

We note that the convex extremals ext{V* {A)) correspond to the pseu- 
dovertices of tropical polytropes [32j (Kleene cones in the min-plus setting), 
and it is known that the number of these may be up to (2(n — 1))!/ {n — 1)! 
pIR [32] . unlike the number of max extremals extQ){V*{A)) which is not 
more than n. 

Max algebraic subeigenvectors give rise to useful diagonal similarity scal- 
ings. A matrix A is called visualized (resp. strictly visualized), if aij = \{A) 
for aU (i, j) G Ec{A), and aij < \{A) for aU (i, j) ^ Ec{A) (resp. a^j < X{A) 
for all {i,j)^E,{A)). 

In the context of max algebra, visualizations have been used to obtain 
better bounds on the convergence of the power method [16], [17] . Strong 
links between diagonal scaling and max algebra were established in . 

Specifically, Corofiary 2.9 of P shows that for a definite A G 
X^^AX is visualized if and only if X = diag(a:) where x is nonnegative 
linear combination of the columns of A* that is positive. 

Strict visualization was treated in a special case [S] , in connection with 
the strong regularity of max-plus matrices. 

A preliminary version of the following theorem appeared in [lOj. 
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Theorem 3.3 Let A G M"^*^ be definite and X = diag(a:) with positive x G 
M" . Then X^AX is strictly visualized if any of the following conditions 
are true: 

1. X is a positive linear combination of all columns of A* ; 

2. A is irreducible andx is a positive log-convex combination of all columns 
of A*. 

Proof. The following argument goes for both cases. In both cases, X 
is positive: for positive linear combinations this is true since a*^ = 1 for 
all and for positive log-convex combinations, Prop. 12.31 assures that A* 
is positive if A is irreducible. As x G ^{A), we have that aijXj < Xi 
for all i^j. By Theorem [23] part 5., for all G Ec{A). If 

{i,j) ^ Ec{A), then, by Theorem [231 part 1., aijZj < Zi for z = A*-, while 
ttijZj < Zi for all z = ^4*^ where /c G [n]. After summing these inequalities 
for all z = A*i^ with positive coefficients, or after raising them in positive 
powers and multiplying, we obtain that aijXj < x^, taken into account 
the strict inequality for z = A*-. Thus x is positive, x^^aijXj = 1 for all 
G Ec{A) and x^^aijXj < 1 for aU {i,j) ^ Ec{A). u 
Note that if A is definite, then every column of A* can be used to obtain 
a visualization of A, which may not be strict. This result was known 
to Afriat P, [3] and Fiedler and Ptak [21], and it has been a source of 



inspiration for many works on scaling problems, see flSl [IHl [23, [33, [SI 
Theorem 13.31 implies the following. 

Proposition 3.4 Let A have X{A) > 0, then there exists X = diag(a:) with 
positive X G such that X^^AX is strictly visualized. 

If A is definite and irreducible then ^4* is irreducible, and in this case 
A* has an essentially unique positive linear algebraic eigenvector, called 
the Perron eigenvector [B]. As it is a positive linear combination of the 
columns of A*, we have the following. 

Proposition 3.5 Let A G M"^" be definite and irreducible and let x be 
the Perron eigenvector of A*. Then X^AX, for X = diag(a:), is strictly 
visualized. 
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We will now give a topological description of strict visualization scalings, 
using the linear hull and relative interior of V*{A). 

By Theorem \2M part 5., for all y G V*{A) and G Ec{A) we have 
ttijijj = Hi. This can be formulated geometrically. For A C consider the 
set 

L{C{A)) = {xeW\ a.jXj = XiA)x, V(z, j) G E,{A)}. 

This is a linear subspace of which contains both V*{A) (as its convex 
subcone) and V{A) (as a max subcone of ^{A)). If B = X^AX with 
X = diag(a:) and x positive, then, by Prop. 12.101 we have C{A) = C{B), 
and we infer that L(C(^)) = {Xy \ y G L{C{B))}. 

Let i^' be a convex cone. The least linear space which contains K will 
be called the linear hull of K and denoted by Lin(i^'). This is a special 
case of the afhne hull of a convex set, see [25] • Denote by By the open ball 
with radius e > and centered at y. The relative interior of denoted 
by n(K), is the set of points y G such that for sufficiently small e we 
have that 5J n Lm{K) C K. If Lin(A') = M", then it is the interior of K, 
denoted by int(i^'). 

The following important "splitting" lemma can be deduced from [12], 
Lemma 2.9. 

Lemma 3.6 Suppose that K C W]_ is a convex cone which is a solution 
set of a finite system of linear inequalities S. Let Si be composed of the 
inequalities of S which are satisfied by all points in K with equality, and 
S2 '•= S\Si be non-empty. 

1. There exists a point in K by which all inequalities in S2 are satisfied 
strictly. 

2. Lin(i^') is the solution set to Si, and n{K) is the cone which consists 
of the points in K by which all inequalities in S2 are satisfied strictly. 

Now we describe all scalings that give rise to strict visualization. 
Theorem 3.7 Let A G and let X{A) > 0. 

1. h{C{A)) is the linear hull of the subeigencone V*{A). 

2. X G n{V*{A)) if and only if, for X = diag{x), the matrix X^^AX is 
strictly visualized. 
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3. i\iy*{A)) contains the eigenvectors of A if and only ifV*{A) = V{A). 

4- If A is definite, then any positive linear combination, and, if A is 
irreducible, also any positive log-convex combination x of all columns 
of A"" belongs to ii{V*{A)) and X^^AX with X = diag(a:) is strictly 
visualized. 

Proof. 1. and 2.: Consider Lemma ESI with K = V*{A), then V*{A) is the 
solution set to the system of inequahties aijXj < x^, and we need to show 
that the inequahties with G Ec{A), and those with (i, j) ^ Ec{A), play 
the role of S*!, and of Lemma [3?6], respectively. For this, we note that by 
Theorem [231 part 6., the inequalities with {i,j) G Ec{A) are satisfied with 
equality for all x G V*{A), and Prop. [3^ implies that there is x G V*{A) 
by which all the inequalities with ^ Ec{A) are satisfied strictly. 

3. : The "if" part is obvious. The "only if" part: from Theorem \1M it 
follows that V*{A) = V{A) if and only if the set of critical nodes is [n]. 
Suppose that ^(^4) is properly contained in V*{A), then there is a node i 
which is not critical. Then for any eigenvector y there is an edge j) for 
which aijyj = yi and obviously {i,j) ^ Ec{A). Hence y ^ n{V*{A)). 

4. : Follows from Theorem [3^ and part 2. ■ 



Note that as V*{A) is the max algebraic column span of A*, its relative 
interior may also contain vectors which are not positive linear combina- 
tions or positive log-convex combinations of the columns of A*. However, 
the relative interior of V*{A), or the set of vectors which lead to strict vi- 
sualization, is exactly the set of vectors that can be represented as positive 
combinations of all cont'ea; extremals in ext(^{V*{A)), see j25] Sect. 2.3. 

We also remark here that a bijection between n{V*{A)) and n{V* (A^)) 
is given by x i-^ since X{A) = X{A^) and if x is positive, then aijXj = 
X{A)xi (resp. aijXj < X{A)xi) holds if and only if aijx^^ = X{A)xJ^ (resp. 
ttijX^^ < X{A)xJ^). In particular, positive linear combinations of rows of 
Kleene stars also lead, after the inversion, to strict visualization scalings. 

If A is strongly definite (that is, X{A) = 1 and an = 1 for all i G [n]), 
then by Prop. 12.61 we have V*{A) = ^(^4), so that ^(^4) is convex and the 
maximum cycle geometric mean can be strictly visualized by eigenvectors 
in n{V{A)). We note that in the case when, in addition, the weights of all 
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non-trivial cycles are strictly less than 1, the strict visualization scalings 
have been described in jHj. 

Strongly definite matrices are related to the assignment problem. By this 
we understand the following task: Given A G M"^" find a permutation tt G 
Pn such that its weight ai^7r(i) • 02,77(2) " • • • " <^n,7r(n) is maximal. A permutation 
TT of maximal weight will be also called a maximal permutation. 

Again, our aim is to precisely identify ("visualize") matrix entries be- 
longing to an optimal solution using matrix scaling. That is, for a matrix 
A with nonzero permutations, find diagonal matrices X and Y such that 
all entries of XAY on maximal permutations are equal to 1 and that all 
other entries are strictly less than 1. 

To do this, we first find a maximal permutation tt and define the corre- 
sponding permutation matrix by 

Using this matrix, we scale A to one of its strongly definite forms {D^)^^A. 
In a strongly definite matrix, any maximal permutation is decomposed into 
critical cycles. Conversely, any critical cycle can be extended to a maximal 
permutation, using the diagonal entries. Therefore, scalings X which visu- 
alize the maximal permutations of {D^)^^A are scalings which visualize the 
critical cycles, and these are given by Theorem 13. 71 After we have done this 
diagonal similarity scaling, we need permutation matrix = (^i7r-i(i)) 
to bring all permutations again to their right place. Thus we get scaling 
x^^{D'^)~^AX which visualizes all maximal permutations. 

Numerically, solving visualization problems by the methods described 
above, relies on the following three standard problems: finding the maximal 
cycle mean, computing the Kleene star of a matrix, and finding a maximal 
permutation. The first problem can be solved by Karp's method [H ESI I2S] , 
the second problem can be solved by the Floyd- Warshall algorithm [31] and 
the third problem can be solved by the Hungarian method [IMj • All of these 
methods are polynomial and require O(n^) operations, which also gives a 
complexity bound for the visualization problems. 

Finally we note that the problem of strict visualization is related to the 
problem of max balancing considered in [ 1351 1371 [35] | . A matrix B is max 
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balanced if and only if each non-zero element lies on a cycle on which 
it is a minimal element. It follows that B is strictly visualized. It was 
shown in [33], EZ), EB] that for each irreducible nonnegative A there is an 
essentially unique diagonal matrix X such that the scaling B = X~^AX is 
max balanced, and hence there is a unique max balanced matrix MB(74) 
diagonally similar to A. Importantly, the matrix MB(74) is canonical for 
diagonal similarity of irreducible nonnegative matrices, that is A is diag- 
onally similar to C if and only if MB (74) = MB(C). A complexity bound 
for max balancing which follows from [351 EH [3B], is 0{n'^), see also [H] 
for a faster version of the max balancing algorithm. 

4 Diagonal similarity scalings which leave a matrix visualized 

Another approach to describing the visualization scalings is to start with 
a visualized matrix and describe all scalings which leave it visualized. 

We first describe the Kleene star of a definite visualized matrix A G 
M"^". Let C*{A) have m strongly connected components C^^, where /i G 
[m], and denote by A^^ the set of nodes in C^^. Denote by Af^,y the (/U-, i^)- 
submatrix of A extracted from the rows with indices in and from the 
columns with indices in N^. Let A^' G W^^'^ be the m x m matrix with 
entries a^j, = maxja^j | i G N^, j G Aj,}, and let E G M"^" be the n x n 
matrix with all entries equal to 1. 

Proposition 4.1 Let A G R"^" be a definite visualized (resp. strictly vi- 
sualized) matrix, let m be the number of strongly connected components of 
C*{A) and let A^ = {a^j_iy), A*j^^ and E^^y be as defined above. Then 

i- ctfifi = 1 for all fi G [m] and a^^ < 1 (resp. a^jy < 1 for fi ^ u), where 

2. any {fi, v)-submatrix of A* is equal to A*^^^ = a*^j^E^jy, where a*^^ is the 
{fi, Tj)-entry of {A'^)* , and E^^y is the {fi, Tj)-submatrix of E. 

Proof. 1.: Immediate from the definitions. 

2.: Take any i G A^, j G A^,, and any path tt = (zi, . . . ,ik) with ii := i 
and ik '.= j. Then tt can be decomposed as tt = ri o di o r2 o . . . o cr^_i o r/, 
where r^, for i G [/], are (possibly trivial) paths which entirely belong 
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to some critical component C^^, with fii := fi and fii = ly, and cr^, for 
i G [/ — 1], are edges between the strongly connected components. Then 
w{7T, A) < w{7t', A), where tt' = t[ o a[ o T2 o . . . o a'i_^ o t[ is also a path 
from i to j such that r/ entirely belong to the same critical components 
as Tj, and (j[ are edges connecting the same critical components as cr^, but 
w{a[^A) = maxjajj | i G A^^^, j G and w{t[^A) = 1. Such a path 

exists, since in a visualized matrix, there exists a path of weight 1 between 
any nodes in the same component of the critical digraph. Thus a*- is the 
greatest weight over all such paths tt'. As tt' bijectively correspond to the 
paths in the weighted digraph associated with the claim follows. ■ 

Note that, after a convenient simultaneous permutation of rows and 
columns, we have that if is a definite visualized matrix, then 





( En 


a^Eu ■ 


■ Oil^Eijj^ 


A* = 


^21^21 


E22 


■ (AnE2m 




y^mlEml 


c4n2Em2 ■ 


Emm J 



(1) 



Note that A^ does not contain critical cycles except for the loops, oth- 
erwise are not the components of C*{A). Hence Ij{A^) = R"^, and we 
can speak of the interior of V*{A^). 

Given a strictly visualized matrix A as above, denote by fi G [m], the 
matrix such that {Ifx)ij = 1 whenever i = j belongs to A^^ and = 

elsewhere, and by A-\- B the direct sum of matrices A and B. 

Theorem 4.2 Let A G R"^"^ be a definite visualized matrix and let m be 
the number of strongly connected components ofC*{A). Let A^ and be 
as defined above. Then X^AX, where X = diag(a:) with x G R" positive, 
is visualized (resp. strictly visualized) if and only if X has the form 

X = Xji + • • • + Xmlm, 

where x is a vector satisfying a^iyXi, < x^ (resp. a^i^x^j < x^), where n ^ v, 
II, v £ [m]. In other words, x G V*{A^) (resp. x G 'mi{y*{A^))). 

Proof. The "if" part: Let X be as described, then the elements a^j, for 
hi € -^/xj do not change after the scaling, so each block A\\ remains 
unchanged, and hence visualized (resp. strictly visualized). For a^j with 
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i G A^^, i G Ny^ /i 7^ z/, we have that aijXj < Xi (resp. cnjXj < Xi), as 

and aij is the maximum over these a^j. Hence X ^AX 
is visuahzed (resp. strictly visuahzed). 

The "only if" part: Suppose that scaling by X leaves A visualized (resp. 
makes A strictly visualized). As A is initially visuahzed, all critical edges 
have weights equal to 1, and x should be such that whenever 
i,j belong to the same A^^^ For i G A"^, j G N^, ^ we should have 
that ). Taking maximum over these a^j, we 

obtain that this is equivalent to a^j^Xj, < x^ (resp. a^j^Xy < x^). 

It remains to apply Lemma 13.61 (with S*! = 0), to obtain that the same 
is equivalent to x G V*{A'^) (resp. x G int(y*(^'^))). ■ 

In the following we discuss some issues concerning linear algebraic prop- 
erties of Kleene cones and Kleene stars. In this context, Kleene stars are 
known as path product matrices, see j2Hl EE I30j . 

For a matrix A G M"^"" with A(^) > 0, we proved that 

L{C{A)) = {x G M" I a,jXj = X{A)x,, G E,{A)}. (2) 

is the linear hull of V*{A). Note that in the case when A is definite and 
strictly visualized, a^j = 1 for all G Ec{A) and X{A) = 1. Also see 
Section 2 for the definition of n{C{A)) and \K{A)\. 

Proposition 4.3 Let A G M';""" have X{A) > 0. 

1. The dimension ofL{C{A)) is equal to the number of strongly connected 
components in C*{A), that is, to n{C{A)) + |A'c(^)|/ 

2. If A IS definite, then C*{A) = C{A*) and h{C{A)) = L(C(^*)). 

Proof. Let A"^, for G [m] where m = n{C{A)) + |A'c(^)|, be the set of 
nodes of C^, a strongly connected component of C*{A). In the case when 
A is definite and strictly visualized, C*{A) = C{A*) is seen from ([T]), where 
< 1 for all n and it is also seen from ([T]) that L{C{A*)) is the 

linear space comprising all vectors x G M" such that whenever i 

and j belong to the same A^^. As L{C{A)) is also equal to that space by 
(EI), we have that L(C{A)) = L{C{A*)). We can take, as a basis of this 
space, the vectors e^, for /i G [m], such that = 1 if j G A^^^ and = 

if J ^ N^, and hence the dimension of L{C{A)) is n{C{A)) + |A^c(-4)|. The 
general case can be obtained using diagonal similarity. ■ 
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Prop. SJ 



enables us to present the following result. 



Theorem 4.4 For any matrix A with X{A) > 0, the max algebraic di- 
mension ofV*{A) is equal to the (linear algebraic) dimension ofL{C{A)), 
which is the linear hull ofV'^{A). 



Proof. It follows from Prop. and Prop. [43] part 1. that both dimen- 
sions are equal to the number of strongly connected components in C{A). 



When A is strongly definite and the weights of all nontrivial cycles are 
strictly less than 1, Theorem 14.41 implies that V*{A) contains n linearly 
independent vectors. This result has been obtained by Butkovic [7], The- 
orem 4.1. One could also conjecture that in this case the columns of A* 
should be linearly independent in the usual sense. However, this is not so 
in general as we show by modifying Example 3.11 in Johnson-Smith [28]. 
Let 

/ 1 5/11 5/11 7/11 7/11 7/ll\ 
5/11 1 5/11 7/11 7/11 7/11 
5/11 5/11 1 7/11 7/11 7/11 
7/11 7/11 7/11 1 5/11 5/11 
7/11 7/11 7/11 5/11 1 5/11 
V7/11 7/11 7/11 5/11 5/11 1 y 

Then the linear algebraic rank of A* is 5, however, by Theorem 14.41 (or 
[Z], Theorem 4.1) the max algebraic dimension of V*{A), and therefore 
the linear algebraic dimension of L{C{A)), are 6. We observe that x = 
[7/11, 7/11, 7/11, 1, 1, 1]^ is a max eigenvector of ^* (hence in V*{A)) but 
it is not in the linear algebraic span of the columns of A*. Finally we note 
that the original form of Example 3.11 in |[2B] provides a Kleene star with 
negative determinant. 



A = A* = 
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